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INTRODUCTION 


There  is  a  growing  realization  that  computer  simulations 
are  valuable  tools  for  investigating  the  characteristics  and 
performance  of  complex  sensor  systems.  In  fact,  many 
simulations  are  used  to  evaluate  proposed  sensor  designs  and 
to  assess  processing  schemes  associated  with  the  sensed  data. 
The  effectiveness  of  such  simulations  depends  on  the  accuracy 
of  the  simulation  models  used  to  represent  the  actual  scene, 
the  sensor  components  and  the  associated  processing.  Fur¬ 
thermore,  the  techniques  used  to  implement  these  models 
affect  the  efficiency  and  costs  associated  with  running  the 
simulation. 

These  considerations  are  especially  significant  in  the 
simulation  of  future  space-based  scanning  sensors  providing 
deep  space  astronomical  measurements(f>2).  These  sensors 
are  expected  to  scan  the  celestial  sky  and  detect  a  number  of 
objects  (e.  g.  ,  stars,  resident  space  objects,  zodiacal  light). 
The  simulation  is  required  to  simulate  the  sensor  detection 
process  with  high  fidelity  and  computational  speed.  These 
conditions  are  difficult  to  satisfy  because  the  simulation  must 
accommodate  high  density  star  fields  present  near  the  galactic 
equator.  In  general,  high  star  density  not  only  increases  the 
amount  of  memory  required  by  the  computer  simulation  but 
also  reduces  the  speed  at  which  simulation  must  operate  to 
achieve  reasonable  accuracy. 

This  paper  presents  a  method  for  simulating  a  multi- 
spectral  scanning  sensor  in  a  fast,  efficient  manner.  The 
models  and  techniques  used  to  accommodate  the  major  radia¬ 
tion  sources  are  discussed,  Ihen  the  specific  structure  and 
capabilities  of  the  simulation  are  presented.  The  computa¬ 
tional  requirements  are  assessed,  and  representative  sample 
outputs  of  the  simulation  are  presented.  The  final  section  is 
devoted  to  conclusions. 

SYSTEM  DESCRIPTION 

The  system  to  be  simulated  is  a  spaced-based,  multi- 
spectral  scanning  sensor  designed  to  collect  information  on 
celestial  objects  and  events.  As  shown  in  Figure  1,  the 
major  on-board  components  of  a  typical  system  include  an 
optics  and  scanner  subsystem,  an  array  of  detector  cells,  a 
set  of  amplifiers  and  filters,  and  a  multiplexer  and  encoder. 
Optical  elements  scan  the  celestial  sky  and  focus  incoming 
radiation  onto  the  array  of  detectors  located  at  the  focal  plane. 
Sensitive  to  a  specific  spectral  region,  the  detector  array  con¬ 
verts  the  radiation  through  the  amplifiers  to  voltage  levels. 

In  turn,  these  amplifier  outputs  are  filtered,  multiplexed  and 
encoded  for  transmission  to  the  ground. 
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Sources  of  radiation  affecting  the  sensor  outputs  include 
point  sources  such  as  stars  and  resident  space  objects,  as 
well  as  extended  sources  such  as  zodiacal  light  and  earth  limb 
radiation.  In  addition,  there  are  impulsive  sources  due  to  the 
incidence  of  nuclear  particles  and/or  radiation  on  the  detec¬ 
tors.  All  of  these  effects  should  be  accommodated  by  the 
simulation. 

SIMULATION  TECHNIQUES 

In  general,  a  simulation  of  this  system  requires  a  de¬ 
scription  of  the  spectral,  spatial  and  temporal  characteristics 
of  the  scene  and  the  sensor.  This  is  exemplified  by  the  power 
incident  on  a  rectangular  detector  by  a  radiation  source, 
which  can  be  expressed  as 


P.j(t)  =  J“  d dx  dy j'f  H.(u.v.X>  G(u-x,  v-y,  X)dudv(l) 
ADj 

where  P^:(t)  denotes  the  power  on  detector  j  by  source  i  at 
time  t,  Hj  denotes  the  spectral  irradiance  at  the  image  plane 
due  to  source  i,  G  represents  the  optical  point  spread  func¬ 
tion,  and  A£)j  represents  the  area  of  dei.ector  j.  The  terms 
x,  y  and  u,  v  denote  the  dependent  variables  representing 
space  coordinates;  X  denotes  the  spectral  wavelength.  Time 
is  implicitly  related  to  the  space  coordinates  by  the  motion 
of  the  scanner.  The  above  integrals  can  be  interchanged 
to  appear  as 

P.  .(t)  dX  f f  H. (u,  v,  X)  dudv  f  f  G(u-x,  v-y,X)dxdy 

U  -too  J  “Leo  ADj  (2) 

The  above  power  expressions  can  be  solved  indirectly  by 
a  Fourier  transform  method^-^.  This  method  involves  solving 
the  convolution  integral  of  (1)  in  the  transformed  spatial  fre¬ 
quency  domain,  relating  the  temporal  frequency  to  the  spatial 
frequency  by  the  scanning  rate,  and  then  inverting  the  re¬ 
sulting  transform  to  the  time  doman,  To  obtain  the  trans¬ 
formed  signal  at  the  filter  output,  the  frequency  transform  of 
the  detector,  amplifier,  and  filter  is  multiplied  by  the  trans¬ 
form  of  Pij(t),  prior  to  the  inversion  process  (assuming 
linear  models  for  these  components).  Then  the  time  domain 
signal  is  obtained  by  inversion.  To  simulate  the  effects  of 
noise  from  background  sources  (e.  g.  ,  zodiacal  light)  and  on¬ 
board  sources,  the  power  spectral  density  (PSD)  of  these 
sources  is  computed  off-line  based  on  the  input  PSD  and  the 
system  transfer  function.  Then  a  random  number  generator 
is  used  to  simulate  noise  corresponding  to  the  second-order 
statistics  of  the  resulting  PSD.  The  availability  of  Fast 
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Fourier  Transform  packages  provides  an  economic  means  to 
implement  this  method(4).  However,  to  simulate  a  large 
array  with  high  resolution  detectors,  the  amount  of  computa¬ 
tion  required  to  achieve  the  desired  accuracy  may  become 
prohibitive. 

A  more  direct  approach  for  implementing  the  above  re¬ 
lations  is  to  make  reasonable  approximations  and  solve  the 
convolution  integral  directly  in  the  space  domain.  Depending 
on  the  amount  of  approximation,  two  basic  methods  can  be 
identified:  a  "brute -force  "  numerical  method,  and  a  semi- 
analytical  method.  Referring  to  (2),  the  numerical  method 
first  evaluates  the  integrals  involving  the  point  spread  function 
G  for  a  given  detector  geometry  during  simulation  initializa¬ 
tion.  This  defines  a  weighting  function  at  every  point  in  a 
two-dimensional  grid.  At  each  time  instant  t,  the  weighting 
function  is  multiplied  by  the  irradiance  Hj  and  the  result  is 
integrated  to  obtain  Pi j (t)  via  a  double  summation  over  an 
area  about  detector  j.  The  final  integration  over  wavelength 
\  is  avoided  by  utilizing  the  average  of  Hi  and  G  over  a  spec¬ 
tral  bandwidth  AX.  When  the  radiation  is  a  point  source,  the 
method  reduces  to  a  simple  multiplication  of  its  intensity  by 
the  weighting  function  to  obtain  the  power  P„(t).  (5) 

The  potential  problem  with  this  numerical  method  is  that 
a  small  grid  spacing  may  be  required  to  achieve  the  desired 
fidelity,  especially  for  simulating  extended  sources.  As  a 
rule,  the  grid  spacing  will  be  some  fraction  of  the  smallest 
detail  of  interest  (e.  g.  ,  the  detector  dimensions,  star  sepa¬ 
ration  distance,  optical  blur  size).  If  the  required  number  of 
grid  points  becomes  large,  the  simulation  may  be  too  slow  or 
expensive  to  run. 

The  semi -analytical  method  avoids  the  above  problem  by 
making  simplifying  assumptions  on  the  radiation  source,  opti¬ 
cal  properties  and/or  array  geometry.  These  assumptions 
permit  the  integrals  in  (1)  or  (2)  to  be  evaluated  analytically. 
As  a  result,  the  incident  power  is  expressed  in  terms  of 
known  functions  which  can  be  efficiently  implemented  in  a 
simulation, 

SIMULATION  ASSUMPTIONS 


The  simulation  under  consideration  uses  the  semi- 
analytical  approach  to  characterize  the  signal  and  noise  of  the 
scanning  system.  The  assumptions  used  to  simplify  the  sensor 
are  given  as  follows: 
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1 .  Gaussian  point  spread  function 


G(x, y, X) 


Dq(\) 


2t r  ct 


-j-e.vp 


±d+x!)| 


2a, 


(3) 


where  a^  denotes  the  optical  blur  size  and  D0  represents  the 
spectral  transmission  function  of  the  optics. 


2.  Bandpass  optical  filter 


Do(\) 


1,  x_ 


<  \  <\ 

2  n  2 


0,  otherwise 


(4) 


where  \n  and  AX  denote  the  nominal  frequency  and  bandwidth, 
respectively. 


3.  Rectangular  aperture  function 


DAj(x.y) 


lx-xCjl  <  W2,  1  y-ycj I  /2 

0,  otherwise 


(5) 


where  w  and  f  denote  the  detector  width  and  length,  respec¬ 
tively;  xcj  and  ycj  define  the  center  location  of  detector  j  in 
the  imageJ  plane. 


4.  Constant  nominal  scan  rate 


x 


•  (t)  = 


CJ 


X  .(t  ) 

Cj  o' 


+  V  t  +  p  (t) 

s  X 


(6) 


y  .(t)  =  y  .(t  )  +  p  (t)  (7) 

cj  'cj  o  y 

where  vg  denotes  the  scan  rate;  px  and  p  define  small 
variations  in  the  detector's  position  due  to  spacecraft  attitude 
motion  or  scanning  jitter;  and  t  denotes  the  initial  reference 
time. 
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5. 


Linear  Detector/Amplifier  Response 


z.  = - —  z.  +  — -  (u.  +  n  .  +  nT.)  (8) 

J  Ta  J  Ta  J  Ej  Ij' 

where  Zj  denotes  the  amplifier  output  of  detector  j;  Ra  and  t& 
represent  the  amplifier  gain  and  time  constant,  respectively; 
Uj  denotes  the  total  signal  incident  on  detector  j  due  to  point 
sources;  njr;j  and  njj  represent  noise  on  detector  j  due  to  ex¬ 
tended  sources  and  impulsive  sources,  respectively. 

6.  Linear  Filter  Response 

q.  =  Aq.  +  bz.  (9) 

where  qj  denotes  an  n  x  1  state  vector  of  the  on-board  filter; 

A  and  b  represent  m  x  n  and  n  x  1  matrices,  respectively. 

SIGNAL  AND  NOISE  RELATIONS 


Sources  of  radiation  are  modeled  in  the  simulation  as 
point  sources,  extended  sources,  or  impulsive  sources.  Point 
sources  which  include  stars  and  resident  space  objects  pro¬ 
duce  the  basic  signals  of  the  system.  Extended  sources  such 
as  zodiacal  light  and  impulsive  sources  such  as  high  energy 
nuclear  particles  are  assumed  to  produce  noise  at  the 
detector. 

The  spectral  irradiance  of  point  sources  can  be  expressed 
at  the  image  plane  as 

IL(u,v,\)  =TloAoHNi(\)  6(u-uoi)  S(v-voi)  (10) 


where 


HNi(M  = 


HSi(\), 

JTi(M 


R 


Star 


Resident  Space  Object 


(11) 


The  terms  T)Q  and  AQ  respectively  denote  the  optical  effi¬ 
ciency  and  optical  entrance  area;  Hgj  is  the  irradiance  of  star 
source  i;  Jq^i  and  R  represent  the  intensity  and  range  of 
resident  space  object  i,  respectively;  6  (•)  denotes  the  Dirac 
Delta  function;  and  uQj  and  vQi  define  the  location  of  the 
point  sources  in  the  image  plane. 
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Combining  Eqs.  (3)  through  (7)  into  (1)  yields  the  incident 
power  on  detector  j  due  to  point  source  i  as 


P..(t)  A  AX  H,,.(X  )  Z(x  .  -  u  .,  w)  Z ( y  .  -  v  I  )  (12) 

i  i  '  ’  A  '  r\  a  1\7 1  '  n  '  '  r  i  ni  '  '  r  1  m 


ij 

where 


'o  o  Ni  n  cj  oi 


cj  oi 


Z(u-v,  a)  =  erf 


fu-v  +  a/2Vcrf  /u-v  -a/2\ 


^  2a . 


\r 2a . 


(13) 


erf(x)  =  -2-  r  exp  (-§2)  d§ 

\I~tt  J 


(14) 


The  signal  Uj  to  the  amplifier  includes  all  incident  power  on 
detector  j  due  to  point  sources,  or 


u. 

3 


£  P«") 


(15) 


where  Np;  represents  the  total  number  of  point  sources 
affecting  detector  j. 

Extended  radiation  sources  are  assumed  to  produce  ran¬ 
dom  photon  noise  at  the  detector.  To  a  first  order  approxima¬ 
tion,  the  statistics  of  photon  noise  depend  on  the  incident 
power  resulting  from  these  sources.  A  measure  of  this  power 
can  be  estimated  by  assuming  that  the  radiance  due  to  these 
sources  is  primarily  unstructured  in  the  scanned  region.  In 
this  case  the  irradiance  Hj  from  extended  source  i  at  the 
image  plane  can  be  expressed  as 

T)  A 

H.(x,y,\)  =— g--  NBi(M 


where  f  denotes  the  focal  length  of  the  optics;  Ngj  represents 
the  unstructured  component  of  the  radiance  from  source  i. 
Substituting  ( 1 6) ,  (3),  and  (4)  into  (1)  yields  the  power  on 
detector  j  due  to  the  extended  source  i  as 


P.. 

i3 


7)  A 

V?  N  . (X  )  A~.  AX 
p  Bi  n'  Dj 


(17) 
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r 


In  this  case,  the  resulting  photon  noise  n  £j  at  detector  j  is 

modeled  as  a  white  Gaussian  random  variable  with  mean  n„ 

y* 

and  auto-correlation  function  R_(t)  given  as  follows: 

E 


n 


E 


Ej 

■£ 

iTi 


P.. 

ij 


(18) 


R£(t)  =  Q£  6  (t-T  ) 

where 


NP) 

2  E 

n 

£ 

P.. 

D 

+£  pij 

i  =  1 
- _ 

i  =  l 

(19) 


(20) 


The  term  En  represents  the  energy  of  a  photon  at  wavelength 
Xn;  T]q  is  the  quantum  efficiency;  and  Njrj  is  the  total  number 
of  extended  sources  affecting  detector  j. 

Impulsive  sources  are  modeled  as  an  impulse  noise  pro¬ 
cess  at  the  detector(6).  In  the  time  interval  (0,T),  the  pro¬ 
cess  consists  of  Njj  independent  pulses  on  detector  j  with  a 
common  waveform,  the  ith  pulse  having  an  amplitude  ay  and 
time  of  occurence  Ty.  Thus,  the  impulse  noise  n..  can  be 
expressed  as 

Si 

nIj  =  £  aijU(t  '  V  •  0StST  (21) 

i=l 


where  u  represents  the  pulse  waveform. 

It  is  assumed  that  the  amplitudes  ay  are  statistically  in¬ 
dependent  random  variables,  characterized  by  a  known  den¬ 
sity  function  f  .  Also,  the  arrival  times  Ty  are  assumed  to  be 
uniformly  and  independently  distributed  in  (0,T).  As  a  result, 
the  number  of  pulses  arriving  in  a  given  time  interval  (0,T) 
follows  a  Poisson  distribution  characterized  by  a  mean  rate  of 
arrival  Xj.  (?)  The  waveform  u  can  be  modeled  by  the  Dirac 
Delta  function  at  the  amplifier  input.  Thus,  given  fa,  X  ,  and 
u,  the  impulse  noise  process  is  specified. 
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To  describe  the  signals  and  noise  for  a  computer  simula¬ 
tion,  a  set  of  difference  equations  is  developed.  Under  the 
above  assumptions,  Eq.  (8)  of  the  amplifier  can  be  integrated 

between  times  t  and  t  ,  ,  to  yield 
n  n+1  ' 

-  A/t 

VW  =  *  a  +  “VW +  wEj(W +  VW  (22> 

where 


wUj(tn+l)  =T~ J.  0 

a  •'t 


t  +A  -(t  +A-t)/t 


m(T)dr 


(23) 


R 


Nii 


V*n-H»=v;  £  aiie 

a  i=l 


-(t  +A-t..)/t 
n  --J  a 


t  .  .  >  t..  (24) 

n+1  ij 


where  A  represents  the  integration  time  ( tn  +  j  - tn ) .  The  term 
wEj  *s  a  w^ite  Gaussian  random  noise  sequence  with  mean  w 
and  variance  a specified  as  * 


w  =  R  n 
E  a  E 


(25) 


4  --  Rl  °E 


[.  -  e-2A^a]  /2r. 


(26) 


In  a  similar  manner,  the  difference  equation  for  the  filter 
can  be  expressed  as 


AA 


qj(tn  +  1>  =  e  <*j<tn)  +  Vj(tn  +  1> 


where 


'A  «>■/" 


t  +A  A(t  +A-t  ) 

n  n  ' 


bz^T  )dT 


(27) 


(28) 


The  integrals  of  (23)  and  (28)  can  be  evaluated  by  making  an 
approximation  on  the  forcing  term  such  as 


9 


(29) 


where 


uj(tn+  1  'T)  SUj(tn+  1}  -  ij(tn+  1)t 
zj(tn+  1  'T)S2j(tn+  1>  -  Zj(tn+  1)t 

S(tn  +  1>  [Uj(tn+  1}  -  uj(tn)] 


Zj(tn+  1>  =i[Zj(tn+  l]  "  Zj(tn>]  (32) 

Substituting  the  above  relations  into  (23)  and  (28)  yields  ex¬ 
pressions  for  wyj  and  vj,  respectively,  which  are  explicitly 
related  to  u:  and  z-  at  times  tn+j  and  tn.  These  difference 
relations  are  readily  incorporated  in  a  digital  computer 
simulation. 

Because  of  the  closed-form  solution  to  the  amplifier  and 
filter  equations,  the  integration  step  size  is  not  limited  by  the 
minimum  time  constant  of  the  system.  Instead,  the  step  size 
is  selected  to  provide  sufficient  accuracy  in  the  approximation 
of  (29)  -  (32).  Reasonable  accuracy  can  be  achieved  by  allow¬ 
ing  a  sufficient  number  of  integration  steps  within  the  detec¬ 
tor's  dwell  time  (w/vg). 

SIMULATION  IMPLEMENTATION 

Based  on  the  models  of  the  previous  sections,  a  digital 
computer  simulation  was  developed  on  the  Control  Data 
Corporation  7600  machine.  As  shown  in  Figure  2,  these 
models  are  grouped  into  an  environmental  block  and  a  sensor 
block.  The  environmental  block  consists  of  a  star  field 
model,  a  zodiacal  light  model,  an  earth  limb  model,  a  resi¬ 
dent  space  object  model,  and  a  nuclear  radiation  model.  This 
block  acts  as  a  driver  for  the  sensor  block,  which  consists  of 
an  optics/scanner  model,  a  focal  plane  detector  model,  an 
amplifier/filter  model  and  a  multiplexer/encoder  model.  In 
addition,  provisions  are  made  to  include  the  effects  of  space¬ 
craft  attitude  jitter,  internal  detector  noise,  and  electronic 
noise. 
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Figure  2.  Block  Diagram  of  a  Scanning  Sensor  Simulation 


The  simulation  can  accommodate  up  to  100  detectors  with 
different  spectral  and  temperal  characteristics.  Also,  pro¬ 
gram  options  exist  to  accommodate  combinations  of  background 
models  and  a  variety  of  amplifiers,  filters,  and  analog  to 
digital  converters.  For  run  efficiency,  logic  is  provided  to 
select  the  set  of  detectors  and  options  of  interest.  Also,  the 
simulation  can  accommodate  multiscanning  operations  involv¬ 
ing  up  to  3000  stars  per  run.  The  memory  required  by  this 
simulation  is  broken  down  as  follows: 


A  rea 

Environmental  Block 
Instruction 
Storage 

Sensor  Block 
Instruction 

Variables,  Parameters 
Storage 

Overhead 

Print,  File  Management 
Plot  Management 


Number  of  Words 


0.  4K 

12.  3K  (Primarily  Star  Data) 


1.5K  (Main)/1K  (Subroutine) 
4.6K 

2.  6K  (Primarily  Cell  Data) 


4.  OK 
4. OK 


The  total  of  30.  4K  words  (60  bit/word)  is  less  than  half  of  the 
small  core  capacity  of  the  CDC  Machine. 


A  major  feature  of  the  simulation  is  that  it  processes  high 
density  star  fields  in  a  fast,  efficient  manner.  First  the  stars 
are  sorted  in  the  order  of  their  arrival  at  the  detectors.  Then 
with  a  known  detector  geometry,  a  set  of  candidate  stars  that 
can  be  detected  at  any  given  time  is  identified.  Furthermore, 
a  masking  technique  is  used  to  allow  only  those  candidate  stars 
within  a  specified  window  about  a  detector  to  contribute  to  its 
output.  By  minimizing  the  number  of  candidate  stars,  the 
sorting  and  masking  operations  reduce  the  amount  of  com¬ 
putation  and  thus  the  run  time.  Finally,  by  implementing 
a  closed  form  solution  of  the  amplifier  and  filter  equations, 
the  simulation  accurately  computes  the  time  history  of  the 
filter  outputs  without  a  large  number  of  integration  steps/ 
run.  In  this  case,  the  maximum  integration  step  size  is 
limited  by  the  detector's  dwell  time  and  not  the  natural 
frequency  of  the  amplifier /filter  model.  By  properly 
sizing  the  step  size  and  the  masking  window,  the  simulation 
provides  a  fast,  accurate  time  history  of  each  detector 
output. 
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Figure  3  shows  the  geometry  of  a  typical  run  involving 
two  linear  arrays  (5  detectors/array)  scanning  a  celestial 
background  consisting  of  stars,  zodiacal  light,  and  earth  limb 
radiation.  The  run  simulates  a  constant  scan  rate  sensor,  a 
band-pass  filter  providing  A-C  coupling,  and  a  bipolar, 
analog-to-digital  converter  providing  logarithmic  compaction 
in  the  encoding  process.  Typical  outputs  of  two  adjacent  de¬ 
tector  channels  are  shown  in  Figures  4  and  5  for  a  field  of 
view  consisting  of  stars  at  a  nominal  density  level  and  five 
times  that  level,  r sspectively.  The  plots  indicate  the  filter 
outputs,  encoder  outputs,  and  the  outputs  of  a  ground  recon¬ 
struction  filter  that  inverts  the  nonlinear  encoding  operation. 

It  is  noted  that  a  high  degree  of  fidelity  exists  between  the  on¬ 
board  filter  outputs  and  the  reconstructed  filter  outputs. 

The  execution  time  for  this  simulation  depends  on  the 
number  of  detector  cells  considered,  the  number  of  point 
sources  scanned,  and  the  number  of  integration  steps/run. 

For  example,  the  execution  time  was  8.  3  seconds  based 
on  10  detector  cells,  100  stars,  and  1254  integration  steps/ 
run.  The  execution  time  increased  to  21  seconds  when  there 
were  500  stars  in  the  field  of  view.  For  the  same  number  of 
stars,  the  execution  time  is  doubled  when  the  length  of  the  run 
is  doubled.  Thus  the  execution  speed  is  more  sensitive  to  the 
total  number  of  integration  steps  than  the  total  number  of  point 
sources  scanned. 

CONCLUSION 

An  approach  for  simulating  scanning  astronomical  sensor 
systems  has  been  presenced.  The  approach  is  based  on  uti¬ 
lizing  simple  but  reasonably  accurate  models  of  the  sensor 
and  its  major  radiation  sources.  Under  reasonable  assump¬ 
tions,  these  models  simplify  the  relations  characterizing  the 
signals  and  noise  in  the  system.  Based  on  these  relations,  a 
highly  flexible  digital  simulation  was  developed.  This  simula¬ 
tion  provides  a  tool  for  assessing  on-board  system  perfor¬ 
mance  and  for  testing  ground  data  processing  software. 

The  simulation  utilizes  sorting  and  masking  techniques  to 
reduce  the  amount  of  computations  associated  with  high  density 
sources.  In  addition,  closed-form  solutions  of  the  differential 
equations  were  implemented  to  provide  high  fidelity  results 
without  a  large  number  of  integration  steps  per  run.  As  a 
result,  the  simulation  is  capable  of  processing  system  outputs 
in  a  fast,  efficient  manner. 
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Figure  4a.  Typical  Outputs  of  Cell  1 
(Nominal  Star  Density) 
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